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^ ■ We present some preliminary results from a series of extremely large, high-resolution N-body 

CN| \ simulations of the formation of early nonlinear structures. We find that the high-z halo mass 

^3 ■ function is inconsistent with the Sheth-Tormen mass function, which tends to over-estimate the 

1^ ' abundance of rare halos. This discrepancy is in rough agreement with previous results based 
on smaller simulations. We also show that the number density of minihaloes is correlated 

l/^ ■ with local matter density, albeit with a significant scatter that increases with redshift, as 

^^ ' minihaloes become increasingly rare. The average correlation is in rough agreement with a 

^^ , simple analytical extended Press- Schechter model, but can differ by up to factor of 2 in some 

' regimes. 



1 Introduction 

The properties of the first luminous objects in the universe remain a big enigma at present 
due to the scarcity of observational data. These objects, the first stars and galaxies, started 
forming very early, only about 100-200 million years after the Big Bang, and their ionizing 
radiation eventually completely reionized the intergalactic medium . This complex process can 
be studied through numerical simulations and semi-analytical models. The former have the 
advantage of being able to describe complex situations and, in particular, the non-linearities of 
the cosmic structures, but are expensive to run and have limited dynamic range. The latter 
are much cheaper to run and thus allow studies of the full parameter space, but inevitably 
involve many approximations and simplifications. The distribution and properties of the low- 
mass halos which host the first stars is an important ingredient in any semi-analytical model 
of reionization. Such small-scale structure potentially has very different properties from larger 
structures we see at later times, as they probe a very different part of the initial power spectrum 
of density fluctuations. Therefore, it is important to check the validity of any models and fits 
to the halo mass function and bias in this new regime. 



2 Simulations 

The results we present in this work are based on series of very large N-body simulations, as 
summarized in Tabled! They follow between 1728^ (5.2 billion) and 5488^ (165 billion) particles 
(the latter at present is the largest cosmic structure formation simulation ever performed) in a 
wide range of box sizes from 2/h Mpc up to 3.2//i Gpc. Spatial resolution ranges from 50 pc/h to 
40 kpc/h, while the particle masses range from 100 Mq up to ^ 6 x 10^^, which yields minimum 
resolved halo masses (with 20 particles) between 2 x lO^M© and lO^^M©. The simulations were 
performed with the code CubeP^MEl*^. 

Table 1: N-body simulation parameters. Background cosmology is based on the WMAP 5-year results. 



boxsize 


Npart 


mesh 


spatial resolution 


"^particle 


^^^halo.min 


2/i-^Mpc 


2048^ 


40963 


48.8/i-^pc 


99.8 Mq 


1996 M0 


6.3/i-iMpc 


1728^ 


34563 


182/i-ipc 


5.19 X 1O3M0 


1.04 X lO^M© 


11.4/i-iMpc 


30723 


6 1143 


186/i-ipc 


5.47 X 103 Mq 


1.10 X lO^M© 


20/i-iMpc 


54883 


109763 


182/i-ipc 


5.19 X 1O3M0 


1.04 X lO^M© 


IM/i-^Mpc 


30723 


6I443 


IMh-^kpc 


5.47 X 10^ Mq 


1.09 X lO^M© 


l/i-^Gpc 


30723 


6I443 


lQ.3h-^kpc 


3.62 X 10^ Mq 


7.24 X 10^° Mq 


3.2/i-iGpc 


40003 


80003 


40.0/i-ikpc 


5.67 X 10^0 Mq 


1.14 X 10^2 Mq 



3 Cosmological structure formation 

In Figure [1] we show the matter (green) and halo (orange) distributions in a thin slice at redshift 
2: = 8 from our largest, 20 Mpc/h box simulation with 5488^ particles. The cosmic web is 
already well-developed and highly nonlinear at these small scales even at such an early time. 
At this point there are over 110 million resolved dark matter halos in the box. The larger, rare 
halos are strongly clustered, with a spatial distribution which is highly biased with respect to 
the underlying density field, and largely follows the high-density filaments and knots. However, 
there are a fair number of smaller halos (minihaloes) which are found in mean and low-density 
(void) regions. The reason for this is that at this time the smallest minihaloes become very 
common haloes, with v — 6c/cr{M) < 1, where 6c ^ 1.69 is the linear overdensity at collapse 
time predicted by the top-hat model and a(M) is the density field variance at the appropriate 
mass scale M. A number of large-scale voids, from a few to ^ 10 Mpc in size, are found in our 
computational volume, as well as a large number of high density peaks. The density therefore 
varies very significantly between sub- volumes. For example, at ^ 0.5 Mpc scale the density 
variation is 1 order of magnitude even at z = 28 and reaches ^2 orders of magnitude at z = 8. 

3. 1 Halo mass function 

We locate the collapsed halos on-the-fly, as the code is running, using a spherical overdensity 
halo finder with overdensity of 178. This is done by first interpolating the particles onto a 
gridded density field (at resolution twice the number of particles per dimension, as listed in 
Table [1]). Local density peaks (with density at least 100 times the average) are located and 
spherical shells are expanded around each peak until the threshold overdensity is reached. The 
resulting object is then marked as a halo (objects with less than 20 particles are discarded as they 
cannot be reliably identified). The halo centre position is calculated more precisely by quadratic 
interpolation within the cell and the particles within the halo virial radius are identified and 



^http://www.cita.utoronto.ca/niediawiki/index.php/CubePM 




Figure 1: Cosmic Web at redshift z = 8 from simulation with boxsize 20 h~^ Mpc and 5488^ = 165 billion particles 
resolving halos with minimum masses of 10^ Mq. Shown are projections of the total density (green) and halos 
(actual size; orange). Slice is 444 h~^kpc thick, images are of the full box (left) and of a zoomed sub-region 

1.8x1.8 h~^Mpc in size (right). 
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Figure 2: Halo multiplicity functions at redshifts z = 30 (left), 
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16.6 (middle) and 2; = 8 (right). 

then the halo properties, e.g. mass, velocity dispersion, center-of-mass, angular momentum, 
radius, etc. are calculated and saved in the halo catalogue. 

The resulting halo multiplicity functions, {M^ /p){dn/dM)^ at 2; = 30, 16.6 and 8.1 are 
shown in Figure [2l Here we conservatively only include well-resolved halos, with at least 50 
(100) particles at z = 8 (higher redshifts). The halo mass functions show significant differences 
from the widely-used Sheth- Tormen (ST) approximation'^ (solid green line), particularly for rare 
halos. At z = 8.1, ST is a reasonably good fit for halos with M^aio < lO^M© (corresponding 
roughly to z/ up to a few) , but over-predicts the abundances of more massive halos by a significant 
factor of up to a few. This is consistent with previous results on the halo mass function at high- 
2;I 4P | Q | ^i higher redshifts the numerical mass functions do not agree with ST by ever larger 
factors, over-predicting the halo abundances by up to an order of magnitude at z = 30. In fact, 
at highest redshift, the classic Press-Schechter mass function (green dotted line) is a better fit, 
although it somewhat underestimates the halo abundances. 

3.2 Halo bias 

The halo mass function is a strongly nonlinear function of the local density. Overdense regions 
behave locally as universe with higher mean density producing exponentially more halos. This 
is directly related to the bias of the halo distribution with respect to the underlying density 
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Figure 3: Number of minihaloes per unit volume vs. local overdensity, 1 + ^, for sub-volumes of size 440/i~^kpc 
at z = 28 (left), z = 20 (middle) and z = 8 (right). Shown are the results based on 2 simulations with same 
resolution, with boxsizes 6.3/i~^Mpc (red points) and 20/i~^Mpc (black points). Also shown are the best mean 
fits for 6.3/i~^Mpc (dark blue line) and 20/i~^Mpc (green line) and the extended Press-Schechter (extPS) model 
predictions (light blue line). Top panel shows the extPS (light blue) and 6.3/i~^Mpc mean (dark blue) in units 

of the respective 20 h~^Mpc box results. 

field and is an important ingredient in many semi-analytical models of structure formation and 
reionization. Such models are also used as sub-grid physics in simulations when they do not have 
sufficient resolution to resolve all halos relevant to the questions being asked. It is therefore very 
important to have a handle on the correlation of mass function with local density. To this end, 
we use data from two of our simulations, with box sizes 20 Mpc//i, 6.3 Mpc/h. Scatter plots 
of the number of minihalos {M^aio < 10^ M©) vs. local density in units of the mean are shown 
in Figure [3l The best-fit mean relations are plotted as well. For comparison we also show an 
analytical bias prescription based on extended Press-Schechter theorj^^. The first observation is 
that results are fairly consistent between the 6.3 Mpc/h and the 20 Mpc/h runs, which have the 
same spatial and mass resolution, but very different volume. Both extremes, high overdensity 
and underdensity, are less well sampled by the smaller box, which is especially evident at high 
redshift, z = 28, but the best mean fits for each box agree with each other reasonably well. 
At later times, z < 20, they become virtually indistinguishable. The analytical model gives a 
relatively good prediction for the correlation at mean density and high redshift, and at high 
density and lower redshift, but can be off by up to a factor of 2 at other regimes. Finally, we 
note that there is a significant scatter in the halo number - local density relation, particularly 
at higher redshift s. At later times the correlation tightens, although mostly in relative terms, 
because the density variation range increases significantly, while the absolute value of the scatter 
remains roughly constant. The origins of this scatter and its effects on semi- analytical and sub- 
grid numerical models are currently under investigation. 
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